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Abstract. We present preliminary results on a new family of distribution functions that are able 
O , to generate axisymmetric, truncated (i.e., finite size) stellar dynamical models characterized by 

toroidal shapes. The relevant distribution functions generalize those that are known to describe 
polytropic spheres, for which all the dynamical and structural properties of the system can be 
00 . expressed in explicit form as elementary functions of the system gravitational potential. The model 

^—l ' construction is then completed by a numerical study of the associated Poisson equation. We note 

that our axisymmetric models can also include the presence of an external gravitational field, such 
7— I ■ as that produced by a massive disk or by a central mass concentration (e.g., a supermassive black 

> 1 hole). 
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CN INTRODUCTION 

cn ! 

Constructing self-consistent collisionless equilibrium models is a key step to understand 
the structure and the dynamics of stellar systems. The construction process usually 
Oh 1 starts with the assignment of a phase-space distribution function (DF) obeying the Jeans 

theorem (and thus a solution of the collisionless Boltzmann equation that describes such 
"gravitational plasma"), which leads to an expression for the density as a function of the 
potential; models are then calculated by solving the Poisson equation, as a non-linear, 
second-order partial differential equation for the potential. In the presence of special 
symmetries (for example, in the case of spherical models) the problem can be reduced 
to the study of an ordinary differential equation, and the model construction is relatively 
straightforward. For more general symmetries the procedure is inherently difficult. In 
general, solutions have to be found numerically and turn out to exist only for a finite 
range of the parameters appearing in the adopted DF 

Here we present the basic properties of the simplest member of a new family of 
DFs that are able to generate axisymmetric, truncated (i.e., finite size) stellar dynamical 
models characterized by toroidal shapes. The complete description of our family of DFs, 
together with the results of N-body simulations aimed at the study of the stability of these 
collisionless equilibrium configurations, will be presented in a separate paper (Ciotti, 
Bertin, Londrillo 2004). We recall that toroidal, self-consistent models are sometimes 
mentioned in the literature (e.g., see Lynden-Bell 1962; the Author expressed concerns 
about their stability). With our models we will investigate this matter further, also 
considering the possible stabilizing effect of a massive dark matter halo. Astrophysical 
applications will address the modeling of peanut-shaped bulges (e.g., see Aronica et 
al. 2003, and references therein) and the problem of the central depression recently 



discovered (Lauer et al. 2002) in the luminosity profiles of some elliptical galaxies 
hosting supermassive black holes. 



THE MODELS 

The particular set of axisymmetric models presented here is defined by the DF 

f(QJz)=fo\J z \ 2a (Q-Qo) p ®(Q-Qo), (1) 

where fo is a (dimensional) physical scale, Q = S — J 2 /{2R 2 ), is the Heaviside step 
function, and a and /3 are two dimensionless constants; the natural coordinates are 
cylindrical, (R,z,q>). In addition, S = - ||v|| 2 /2 = ^x - (v 2 R + v\ + v 2 )/2 is the 
binding energy per unit mass, J z = Rvy is the axial component of the angular momentum 
(per unit mass), and Qq and 7? a are a truncation energy and an anisotropy radius, 
respectively. The function = ^tC^z) represents the total gravitational potential; 
for the moment, we only require that ~ 0(r _1 ) for r — > °°, where r = \/R 2 +z 2 . 
We allow for an external component also, i.e., = l I / + l I / ext> where *P is the potential 
associated with the density distribution derived from eq. (1) as given by the following eq. 
(4), while ^ext is taken to be a given function (equal to zero in the fully self-consistent 
case). Note that we can write 

,,2 v 2 / r,2 
v m v v I ^ R 



e ^ T - f - f[l + -j, (2) 

where = v\ + v\ is the square of the meridional (or poloidal) velocity component. 
Some dynamical properties of the systems associated with eq. (1) are the following: 
(/) no net rotation is present, and (ii) the two components o\ and c z 2 of the velocity 
dispersion tensor are equal. In addition, note that for a = if we formally consider 
R a — > oo we recover the polytropic spheres, while for integer a and /3 the distribution 
function is a sum of Fricke's (1952) DFs. The functional dependence of the system 
structural and dynamical properties on ^x is obtained by integration over £2(x), the 
velocity section of phase-space over which / > at given x. The coordinates (v, £, to 
be used to integrate the various quantities over the rotationally symmetric (ellipsoidal) 
region Q(x) are naturally defined by 

• s- f ■ c ■ f vcosC 

v R = vsrn£cos£, v z = vsinCsin<^, V<p = —====, (3) 

where 0<£<7T, < ^ < 2k. With this choice, the ellipsoidal radius v of the velocity 
section Q(x) runs in the range < v < a/2[ v I / x(x) — Qq]. Integration over the two angular 
coordinates gives 

p = 4Kf 0g (R) B{a + 3 A ^ + l )2g+1/2 CP T ~ eo) a+ ^ +3/2 e(^T - Qo), (4) 



where B{x, y) = j$t x 1 ( 1 — t) y l dt is the Euler Beta function, 



R 2a 

8{R) ~ (l+/?2//?l)« +1 / 2 ' 



and for convergence it is required that a > — 1/2, j8>— 1. Moreover, 

- ^w ^^^r ^ - eor*-/^ - eo) 



(2a + l)(2a + 3) 

and 



(6) 



K- 5;y (a+5/ ^' )yt!,i ^- g »» ww; ^- ft »' m 

The number of free parameters (excluding those associated with the external potential) 
needed to describe a model completely is five: the three dimensional quantities /o, 
i? a , and Qq, and the two dimensionless constants a and /3. The relation between the 
meridional velocity dispersion and the azimuthal (toroidal) velocity dispersion o^, 
can be described by the anisotropy distribution 

* Z) = 1 "^- R 2 /R 2 + l ■ (8) 

The velocity dispersion tensor is isotropic when a = 0, tangentially anisotropic when 
a < 0, and meridionally anisotropic when a > 0. Note that the velocity dispersion 
anisotropy is constant on cylinders. It can be proved that eq. (8) holds also when the 
factor (<2 — Qo)P in eq. (1) is replaced by a generic function h(Q — Qq) (Ciotti et al. 
2004). We now impose the system consistency, i.e., we require that 

A^=-47tGp, (9) 

with ¥ a nowhere negative function. In order to obtain the numerical solution of eq. (9), 
the Poisson equation is first recast in dimensionless form, by referring to the physical 
scales R a for lengths and Qq for potentials. After rescaling, we are then led to the solution 
of the following problem: 

A0 = -Ap, p = g(R)(<fc-iy®(<fr-l), (10) 

where y=a + /3 + 3/2>0, and the dimensionless quantities are defined as 0t = / <2o> 
R =R/R a , and g(R) = R 2a /(l + R 2 ) a+1 / 2 , with 

. _ l6n 2 GfoRl a + 2 Q^ +1/2 B(a + 3/2^ + l)2 a +^ 2 

x = 2^Ti • 

It follows that 

p = -^£ (12) 
P GR 2 An K ] 




FIGURE 1. Isodensity contours of a meridional section of the dimensionless density distribution asso- 
ciated with the models (a, 7) = (1,3) for X ~ 3.1 (left) and X ~ 10.3 (right). 



«-™<*Z™rz!> (1 3) 

Equation (10), a classical non-linear elliptic partial differential equation, is then solved 
by using a Newton iteration scheme, described in detail by Ciotti et al. (2004). In 
general, for assigned model parameters we found convergence only for finite intervals 
of A. In Figure 1 we show the meridional sections of the density distribution of two 
representative truncated and fully self-consistent (^ext = 0) toroidal models, obtained 
by fixing a > (see eqs. [5] and [10]), for two different values of A. 
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